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I .   INTRODUCTION 
1.1  General  Remarks 

The  purpose  of  the  research  presented  here  was  to 
develop  an  automatic  detection  scheme  for  signals  resembling 
those  encountered  in  Raman  spectroscopy.  The  detection  of 
Raman  spectral  peaks  presents  a  special  problem  because  of 
the  low  signal-to-noise  ratios  (SNRs)  that  are  often  dealt 
with  and  because  the  shapes  of  the  peaks,  which  are 
Lorentzian  or  Gaussian  in  form,  make  them  difficult  to 
detect. 

Various  methods  have  been  used  for  peak  enhancement. 
Matched  filtering,  ensemble  averaging,  and  boxcar  averaging 
are  different  methods  for  improving  signal-to-noise  ratios 
of  signals.  A  matched  filter  is  a  filter  whose  impulse 
response  is  a  time-reversed  version  of  the  profile  of  the 
expected  signal.  A  matched  filter  can  be  shown  to  maximize 
the  signal-to-noise  ratio  when  the  noise  is  white.  However, 
the  effectiveness  of  the  filter  depends  on  how  closely  the 
signal  matches  the  filter  response.  Ensemble  averaging 
consists  of  making  a  number  of  runs  and  averaging  the 
resulting  sequences .  When  N  runs  are  used  the  SNR  is 
improved  by  a  factor  of  Vn".  The  problem  with  ensemble 
averaging  is  that  a  longer  time  is  needed  to  obtain  results. 
In  boxcar  averaging  the  data  sequence  is  divided  into  groups 


of  n  adjacent  points  and  the  points  in  each  group  are 
averaged  The  SNR  ratio  is  improved  by  a  factor  of  Vn". 
However,  an  improvement  in  the  SNR  requires  decreased 
resolution. 

The  adaptive  peak  detector  [1]  uses  an  adaptive  filter 
to  detect  changes  in  the  statistics  of  the  input  which  would 
indicate  the  presence  of  a  signal.  It  has  been  demonstrated 
to  be  capable  of  detecting  peaks  at  SNRs  of  from  1.0  down  to 
0.3;  the  main  problem  is  that  there  are  several  parameters 
that  must  be  adjusted. 

1.2  Summary  of  Principal  Results 

The  general  strategy  was  to  first  enhance  the  peaks  by 
use  of  a  filter  matched  to  a  typical  pulse  shape,  and  to 
then  follow  this  with  an  automatic  detector.  The  detector 
was  developed  by  an  application  of  statistical  decision 
theory.  The  decision  criteria  chosen  for  study  were  the 
Neyman-Pearson  test  and  the  Wald  (or  sequential)  test, 
mainly  because  they  do  not  require  knowledge  of  the  a  priori 
probabilities  of  the  signal.  The  Wald  test,  which  has  the 
feature  of  using  a  variable  number  of  samples  to  make 
decisions,  was  chosen  for  implementation. 

An  addition  to  the  scheme  that  was  found  to  be 
necessary  was  to  fit  and  remove  a  trend  from  the  output  of 
the  matched  filter,  so  that  the  threshold  comparison 
procedure  in   the  detector  could  work  properly  in  identi- 


fying  individual  peaks  that  are  closely  grouped. 

The  detection  scheme  has  some  limitations  in  that  an 
estimate  of  the  noise  variance  must  be  available,  and  the 
performance  of  the  filter  depends  on  how  well  the  peaks 
match  its  shape.  Some  parameters  involved  in  the  sequential 
detector,  particularly  the  theoretical  probabilities  of 
false  detection  and  missed  detection,  have  an  effect  on  the 
detector  output.  However,  it  may  be  possible  for  the 
detector  to  work  without  need  for  adjustment  of  these 
parameters  on  the  part  of  the  observer. 

1.3  Structure  of  the  Thesis 

The  theoretical  background  for  the  detector  is  given  in 
Section  II.  This  begins  with  a  general  discussion  of 
hypothesis  testing,  followed  by  a  detailed  presentation  of 
the  Neyman-Pearson  and  Wald  (sequential)  tests,  including 
the  development  of  these  tests  for  the  case  of  an  unknown 
signal  in  Gaussian  noise,  the  situation  of  interest  here. 

Section  III  presents  some  details  of  how  the  detection 
procedure  and  the  testing  of  it  are  performed  in  software. 
Included  are  a  discussion  of  the  generation  of  the  data 
sequences,   the  matched  filter,  and  the  detector. 

Section  IV  discusses  first  how  tests  on  several  data 
sequences  were  used  to  determine  how  the  performance  of  the 
detector   depends   on  the  selection  of   the   detector 


parameters.  Then  a  set  of  tests  is  presented  and  some 
conclusions  are  drawn  about  the  detector's  performance  in 
terms  of  its  ability  to  detect  pulses  at  low  SNRs  and  its 
ability  to  avoid  false  detections. 


II.   MATHEMATICAL  BACKGROUND 
2.1  Statistical  Decision  Theory 

Statistical  decision  theory,  or  hypothesis  testing,  was 
first  developed  by  Neyman  and  Pearson  and  later  expanded  by 
Wald  and  others  [2].  The  purpose  of  hypothesis  testing  is 
to  test  two  or  more  alternative  hypotheses  and  to  determine 
which  is  the  best  choice,  according  to  some  criterion  of 
optiraality.  Some  of  the  more  common  examples  of  decision 
criteria  include  the  Bayes,  the  Ideal  Observer,  and  the 
Neyraan-Pearson  criteria.  What  follows  is  a  discussion  of 
decision  theory  in  terms  of  its  application  to  the  problem 
of  signal  detection. 

In  a  detection  problem  the  first  hypothesis  (denoted 
Hq)  is  that  noise  alone  is  present,  and  the  second 
hypothesis  (H, )  is  that  signal  plus  noise  is  present.  The 
situation  discussed  here  is  actually  an  example  of  composite 
hypothesis  testing  [ 3 ] ,  where  the  second  hypothesis  is 
composite;  i.e.,  one  or  more  parameters  of  the  signal  (such 
as  the  amplitude)  may  have  a  range  of  possible  values  rather 
than  a  single  value.  There  may  be  more  than  two  hypotheses 
if  a  decision  is  to  be  made  from  among  various  classes  of 
signals;  a  binary  detection  problem  is  one  in  which  there 
are  two  hypotheses. 

The  set  of  signals  that  are  possible  is  called  the 
signal  space  (Fig.   2.1),  while  the  observations  that  are 


made  are  considered  to  be  elements  of  the  observation  space. 
An  observation  will  consist  of  one  or  more  samples  of  the 
received  signal,   and  when  there  is  more  than  one  sample  the 

observation  may  be  represented  as  a  vector  z  =  [z,   z-   ... 

T 
z  ]   of  the  individual  samples.   For  a  binary  test  the 

observation  space  is  divided  into  two  decision  regions  Z^ 

and  Z^,   corresponding  to  the  two  hypotheses.    Some  kind  of 

decision  criterion  must  be  used  to  decide  which  hypothesis 

should  be  chosen  for  any  given  value  of  z,   that  is,   to 

decide  the  boundary  between  the  decision  regions. 


Decision  Di 


Decision  Do 


Figure  2.1.   Signal  space  and  observation  space. 

In  some  cases  application  of  one  of  these  criteria 
leads  to  a  likelihood  ratio  test.  The  likelihood  ratio7V(z) 
is  a  ratio  of  conditional  probabilities,  i.e. 

A(z)  =  p(z|H^)/p(ziHQ)  (2-1) 

A  likelihood  ratio  test  consists  of  comparing  this  ratio  to 
some  threshold  value  so  as  to  decide  from  among  the 


hypotheses.  A  larger  liklihood  ratio  means  that  H,  is  more 
likely  to  be  true. 

In  a  binary  detection  problem,  there  are  two  possible 
kinds  of  errors:  deciding  'signal  present'  when  there  is  no 
signal  (Type  I  error)  and  deciding  'signal  absent'  when  a 
signal  is  present  (Type  II  error). 

Some  tests  require  knowledge  of  the  a  priori  signal 
statistics,  and  may  also  require  that  costs  be  assigned  to 
each  type  of  error  and  to  each  type  of  correct  decision 
(signal  present  or  signal  absent).  For  example,  the  Bayes 
criterion  requires  an  assignment  of  all  costs,  and  the 
Minimiom-Probability-of -Error  criterion  is  a  special  case  of 
the  Bayes  criterion  in  which  the  cost  of  either  type  of 
error  is  1  and  the  cost  of  a  correct  decision  is  0  [4].  Two 
tests  which  are  useful  when  this  information  is  not 
available  are  the  Neyman- Pear son  test  and  the  Wald  test. 

2 . 2  Neyman-Pearson  Test 

The  Neyman-Pearson  test  does  not  require  knowledge  of 

the  a  priori  probabilities  because  it  is  only  concerned  with 

the   conditional  probabilities  of  error   p(D^JHq)   and 

P(Dq|H-),  and  not  the  total  probabilities.   The  objective  of 

the  test  is  to  keep  both  the  probability  P,  of   false 

detection  and  the  probability  P  of  a  miss  as  low  as 

m 

possible.   Since  decreasing  one  of  these  probabilities  is 


likely  to  mean  increasing  the  other,   the  strategy  that  is 

used  is  to  select  a  value  for  P^  and  to  then  minimize  P  . 

t  m 

This  is  the  same  as  maximizing  the  probability  P,  of 
detection  where  P^  =  1  -  P^^.  The  miss  probability  is 
minimized  subject  to  the  constraint  P^  =  a',  where  a'  is  the 
chosen  value  for  P^,  using  the  calculus  of  extrema  [5].  From 
this  is  developed  a  likelihood  ratio  test  where  the  decision 
strategy  is: 

P(z|H^)/p(zjHq)  <  A    decide  Hq  (2-2a) 

P(z|H^)/p(z|Hq)  >  7\    decide  H^  (2-2b) 

where  7\  is  the  threshold  of  the  test. 

The  threshold  is  determined  by  the  choice  of  a'  and  is 
calculated  from 

a'  =  /p(zjHQ)dz  (2-3) 

The  noise  process  in  Raman  spectroscopy  can  be 
considered  to  be  Gaussian  with  a  mean  of  zero.  If  a  known 
variance  j^  is  assumed,  then  for  a  single  sample  z  the 
conditional  probability  distribution  functions  for  the  two 
hypotheses  will  be: 

P(z|Hq)  =  (l/V^77<r)exp{-z^/2<y^}  (2-4a) 

p(z|H^)  =  (l/V^nor)exp{-(z-m)^/20^}  (2-4b) 

where  m  is  the  signal  amplitude. 

The  likelihood  ratio  then  becomes: 

A(z)  =  exp{(-(z-m)^  +  z^)/2a^} 

=  exp{(2mz-m^)/2ff^}  (2-5) 


It  is  conimon  to  work  with  the  natural  logarithm  of  the 
likelihood  ratio,  L(z),  so  that  the  test  becomes: 

L(z)  =  (2mz-m^)/20^  <  ln(7^)     decide  Hq  (2-6a) 

L(z)  >  ln(A)  decide  H^  (2-6b) 
or 

z  <  [2<J^lnCX)  +  m^]/2m  =  z^      decide  Hq  (2-7a) 

z  >  z^                        decide  H^  (2-7b) 
The  threshold  z^  is  determined  from  the  probability  of 
false  detection  as  shown  in  Egn.  (2-3): 

M 

a'  = /(l/\/27r(r)exp{-z^/2or^}  dz  =  erfc(z./(T)  (2-8) 
Thus,  for  a  chosen  value  of  a'  the  threshold  can  be 
calculated  from  the  inverse  erfc  function.  Note  that  since 
z^  =  cr-erfc   (a'),  the  noise  variance  must  be  known. 

2.3  Neyman-Pearson  Test  With  Multiple-Sample  Observations 

Usually,  a  decision  is  based  on  a  set  of  samples  rather 
than  a  single  sample.  In  this  case  we  work  with  the  joint 
conditional  pdf's,   e.g.   p(z^,  Z2, . . • ,z^ JHq)  for  k  samples. 

The  set  of  samples  can  be  treated  as  a  vector  z  =  [z^  z^  ... 

T 
Zj^]  .   The  observation  space  is  then  a  k-dimensional  vector 

space  which  is  divided  into  decision  regions  Z^   and  Z,  by  a 

decision  surface. 

Again,   for  the  case  of  an  unknown  signal  in  Gaussian 

noise,  the  joint  conditional  probabilities  would  be 

P(zjHq)  =  l/(27Ta2)^/2  exp{-z'^z/2cJ^}  (2-9a) 

p(z|Hj^)  =  l/(27Tff2)^/^  exp{-(z-m)'^(z-m)/20^}     (2-9b) 


T 
where  ra  =  [m^  m2  ...  mj^]  ,  and  the  likelihood  ratio  is 

jV(z)  =  exp{{m'^z  -  .5m'^m)/ar^}  (2-10) 

The  decision  rule  is  then 

m^z  >  ff^lnOv)  +  O.Sm'^m     decide  H^  (2-lla) 

m'^z  <  <r^ln(>k)  +  O.Sm'^m     decide  Hq  (2-llb) 

It  can  be  seen  that  some  knowledge  of  m,  i.e.  of  the  shape 
of  the  signal,  may  be  required,  or  at  least  some  estimate 
must  be  made.  This  will  also  be  seen  when  the  Wald  test  is 
developed  later. 

If  a  constant  but  unknown  signal  amplitude  is  assumed, 
then  a  test  is  developed  comparable  to  that  of  Eqns.  (2-7) 
and  (2-8),  involving  the  sum  of  the  sample  values.  This 
result  will  be  referred  to  in  Section  III. 

It  may  be  difficult  to  decide  what  sample  size  will  be 
best  when  the  signals  are  not  precisely  known.  In  such 
cases  it  might  be  more  useful  to  employ  a  modified  form  of 
the  Neyraan-Pearson  test  known  as  sequential  testing,  in 
which  the  sample  size  varies  according  to  the  number  of 
samples  required  to  reach  a  decision. 

2.4  Sequential  Detection   (Wald's  Test) 

In  the  sequential  detection  procedure  the  nimiber  of 
samples  per  decision  is  not  fixed,  but  varies  according  to 
number  needed  to  reach  a  decision  with  the  desired  probabi- 
lities of  Type  I  and  Type  II  errors.   The  test  involves 
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comparing  the  joint  likelihood  ratio  of  all  the  samples 
taken  with  two  thresholds.  If  the  likelihood  ratio  is  above 
the  higher  threshold,  a  signal  is  declared  present;  if  it  is 
below  the  lower  threshold,  it  is  decided  that  no  signal  is 
present;  and  if  the  likelihood  ratio  is  between  the  two,  no 
decision  is  made  and  another  sample  is  taken. 

If  the  samples  are  statistically  independent,   then  the 

joint   conditional   density  function  p(z.|H,)    can   be 

i  ■■        J 

calculated  as  FiP^^il^v^  >    and  therefore  7\.( z  . )  =  £1(2.). 

The  desired  probabilities  for  both  Type  I  and  Type  II 

error  must  be  specified.    If  we  assign  P^  =  a'  and  P  =  b' 

r  m 

then  the  thresholds  for  the  likelihood-ratio  test  that  were 
derived  by  Wald  [6]  are 

n^  <  (l-b)'/a'  (2-12a) 

and 

Hq  >  b'/(l-a')  (2-12b) 

where  n^  is  the  upper  threshold  and  nQ  is  the  lower 
threshold. 

The  inequality  relationships  arise  from  the  require- 
ments that  A(Zj)  >  n^  andA(z.)  _<  n^  for  acceptance  or 
dismissal,  respectively,  of  a  signal  [7].  If  it  can  be 
assumed  that  the  log-likelihood  ratio  L(2.)  has  small 
increments  L(z.),  then  the  thresholds  can  be  defined  by  the 
equalities 

In(n^)  -  ln[(l-b')/a']  (2-13a) 
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and 


In(nQ)  -   ln[b'/(l-a' )]  (2-13b) 


The  average  number  of  samples  required  to  reach  a 
decision  under  each  hypothesis  has  been  calculated  [8],  and 
it  has  been  shown  that  Wald's  test  will  on  the  average 
require  fewer  samples  to  detect  a  signal  than  a  fixed-sample 
test  with  the  same  probabilities  of  error. 

The  form  of  the  Wald  test  for  the  case  of  an  unknown 
signal  in  Gaussian  noise  will  now  be  developed.  The  noise 
pdf  is  again  described  by  a  mean  of  zero  and  a  standard 
deviation  of  <T.  The  density  functions  for  a  single  sample 
z.  are  as  in  Eqns.  (2-4),  with  a  mean  of  m. ,  and  the  joint 
conditional  density  functions  for  the  jth  sample  are 

p(z.iHn)  =  [l/(2rTor^)^/2]exp{-(l/2<y^  t(z.)^)}     (2-14a) 


2 


and 

P(z.|Ht)  =  [l/(2Trci^)^/2]exp{-(l/2or^)  t'^z.-m.] 
■J   ■"■  (.1  ■■-   -■■ 

(2-14b) 
and  the  joint  likelihood  ratio  is 

A(z.)  =  exp{-(l/2ff^)  t[(z.-m. )^  -  (z.)^]} 

=  exp{-(l/2<r^)  i;[2m.z.  -  (m.)^]}  (2-15) 

the  likelihood  ratio  test  could  then  be  expressed  as 

ym.z.  <  cr^ln(n„)  +  0.5  >*(m.)^     decide  H^     (2-16a) 

2Jm.z.    >  a^ln(nT)  +  0.5  ^"{m.)^     decide  H,      {2-16b) 
;»i  -'■  -1-  -L         ill   1  -L 

If  it  is  assumed  that  all  the  m.  are  identical,  i.e.,  the 

12 


signal  amplitude  is  constant,  then  this  can  be  simplified  to 
Yz^  <  (0'^/m)ln(nQ)  +  0.5jm  decide  Hq  (2-17a) 
fz.    >    (<r^/m)ln{n^)    +   0.5jm     decide  H^        (2-17b) 

with  the  sum  of  the  samples  being  compared  to  the  thresholds 

on  the  right-hand  sides  of  Eqns.  (2-17). 
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III.   SIMULATION  OF  THE  DETECTION  PROCEDURE 
A  program  was  written  to  simulate  the  matched  filter 
and  the  detector  and  to  generate  signals  for  testing.   The 
complete  program  listings  are  in  the  Appendix. 

Signal  Generation.  Two  forms  for  the  signal  pulses 
that  are  of  interest  in  Raman  spectroscopy  are  the 
Lorentzian  and  the  Gaussian  forms,  since  the  Lorentzian 
shape  is  typical  of  the  spectra  of  solids,  and  the  Gaussian 
shape,  for  liquids.   The  Gaussian  form  is 

s(V)  =  A  exp{-(V-VQ)2/4in(2)«^}  (3-1) 

where  v  is  the  wavenumber,  Vq  is  the  wavenumber  of  the 
center  of  the  peak,  and  a  is  the  half -width  at  half -height 
(hwhh)  of  the  pulse.   The  Lorentzian  form  is: 

s(V)  =  A/[l+(V-VQ)^/flc2]  (3-2) 

For  the  Gaussian  shape  the  relationship  of  the  hwhh  to 
standard  deviation  C  is  C  =  oc/V21n2. 

The  program  allows  the  generation  of  two  types  of 
signal  sequences.  In  one  type,  the  heights,  widths,  and 
spacings  of  the  pulses  are  all  randomly  chosen,  being 
uniformly  distributed  between  selected  maximum  and  minimum 
values.  For  position  and  hwhh  these  are  set  within  the 
program.  The  maximum  amplitude  is  specified  by  the  user,  and 
the  minimum  amplitude  is  one-half  of  the  maximum.  The 
limits  now  used  are  50  to  150  wavenumbers  (or  cm"-*")  for  the 
spacing  between  the  centers  of  pulses,    and  5   to  20 
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wavenumbers  for  the  hwhh.  Examples  can  be  seen  in  Figs.  4.9 
and  4.17.  In  the  second  type  of  sequence,  the  amplitude  and 
spacing  are  constant,  although  the  width  varies,  advancing 
in  even  steps  from  5  to  20  wavenumbers  for  the  set  of  10 
pulses  (Fig.  4.1) . 

After  both  the  signal  and  noise  have  been  generated, 
the  signal-to-noise  ratio  (SNR)  for  each  pulse  is 
calculated.   The  definition  of  SNR  used  here  is 

SNR  =  ^s(Vi)'^/i:n{V/)  .  (4-3) 

The  estimate  for  the  variance  of  the  noise  is  the  mean- 
square  noise  level,  calculated  as 

ff^  =  (1/N)  Xn{Vi)2  (4-4) 

Matched  Filtering.  Before  the  statistical  decision 
procedure  is  carried  out,  the  input  signal  is  passed  through 
a  filter  that  is  matched  to  a  Lorentzian  pulse  of  typical 
hwhh. 

The  discrete-time  impulse  response  of  the  filter  is  128 
points  (wavenumbers)  in  length  with  the  pulse  peak  at  63, 
and  it  has  a  maximum  amplitude  of  1.0  and  a  hwhh  of  10.0. 

The  assumption  is  made  that  the  noise  component  of  the 
signal  at  the  filter  output  has  a  Gaussian  density  function. 
This  assumption  should  hold  if  the  filter  is  linear.  It  has 
been  verified  that  the  noise  at  the  filter  output  is  roughly 
Gaussian  by  computing  histograms  of  the  output  of  the  filter 
resulting  from  an  input  of  noise  alone. 
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The  process  of  filtering  involves  performing  a  discrete 
convolution  of  the  filter  response  and  the  signal.  The 
overlap-add  method  of  fast  convolution  is  used  [9]. 

High-Pass  Filtering  (Trend  Removal).  A  problem  that 
appeared  in  early  tests  was  that  if  pulses  are  so  close  that 
in  the  matched  filter  output  the  signal  level  does  not 
decrease  far  enough,  the  detector  may  not  resolve  them  as 
separate  pulses.  One  idea  for  overcoming  this  problem  is  to 
filter  the  output  so  as  to  remove  very  low-frequency 
components,  since  large  pulses  that  are  close  together 
create  a  general  trend  which  rises  during  the  occurrence  of 
such  pulses. 

This  was  done  by  having  the  program  write  the  output 
sequence  from  the  matched  filter  to  a  file,  and  then  using 
RALPH,  a  general-purpose  signal-processing  program  [10],  to 
remove  a  polynomial  trend.  A  fifth-degree  polynomial  was 
used.  The  effect  can  be  seen  by  comparing  Figs.  4.3  and  4.4 
or  Figs.  4.11  and  4.12.  The  result  was  that  the  large  peaks 
were  lowered  enough  so  that  they  could  be  separated  more 
easily.  One  problem,  as  will  be  seen  in  the  results,  is  that 
smaller  pulses  are  sometimes  lowered  so  much  that  they  are 
more  difficult  to  detect. 

Sequential  Detector.  The  algorithm  for  carrying  the 
sequential  detection  procedure  is  shown  in  Fig.  3.1. 

The  program  originally  allowed  the  sampling  to  continue 
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S  =  Sum  of  Samples    t.  Zi, 

n'o  =  Lower  Threshold  of  Wald   Test 

n;  =  Upper  Threshold  of  Wald  Test 

T  =  Threshold  of  Neyman- Pearson  Test 

N  =  Maximum  Sample  Length 
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Figure  3.1.   Sequential  detection  algorithm. 

indefinitely  until  a  decision  was  reached.  However,  this 
might  cause  a  problem  in  that  the  thresholds  are  increased 
for  each  additional  sample,  and  a  peak  that  appears  after 
many  samples  have  been  taken  may  be  harder  to  detect.  A 
good  suggestion  for  a  maximum  sample  length  might  be  the 
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half-width  at  half -height  or  the  width  at  half -height  of  a 
typical  pulse,  and  a  maximvun  length  of  15  samples  is  now 
used.  Once  this  maximum  length  is  reached,  a  decision  would 
be  made  on  the  basis  of  the  samples  already  taken. 

This  decision  is  made  using  a  multiple-sample  Neyman- 

Pearson  test,   with  a  probability  of  false  detection  of  a'  = 

-4 
10   .   When  a  constant  but  unknown  signal  amplitude  is 

assumed,  the  test  shown  in  Eqns.  (2-11)  can  be  developed  so 
that  the  sura  of  k  samples  is  compared  to  the  threshold  T^  = 
k  erfc"  (a')o-[ll].  For  k  =  15  and  a'  =  lo"^  this  becomes  T^ 
=  3.8V15  cr  . 
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IV.   EXPERIMENTAL  RESULTS 
4.1  Effect  of  Detector  Parameters 

The  sequential  detection  procedure  presented  above 
involves  a  requirement  that  some  parameters  be  selected.  In 
setting  the  threshold,  the  theoretical  probabilities  of 
error  must  be  decided  upon,  as  must  an  estimate  for  the 
signal  amplitude,  as  shown  by  Eqns.  (2-12)  and  (2-17).  The 
rate  of  sampling  must  also  be  chosen.  It  is  desired  that 
the  detector  should  give  consistent  results  on  pulses  that 
vary  in  amplitude  or  width  when  a  single  set  of  parameters 
is  chosen.  Some  discussion  of  these  parameters  follows. 

Signal  Level.  The  sequential  test,  as  described  in 
Section  II,  is  based  on  the  assumption  of  a  known  signal, 
and  an  assximed  signal  level  is  used  in  setting  the  threshold 
as  shown  in  Eqns.  (2-17).  Any  signal  which  exceeds  this 
level  will  be  detected  sooner  and  a  signal  which  stays  below 
this  level  will  not  be  detected  as  soon  and  will  have  a 
lower  probability  of  being  detected.  In  the  tests  it 
appeared  that  a  higher  assiomed  signal  level  actually  seemed 
to  help  in  detection  of  weak  signals,  and  also  made  false 
detections  more  likely  in  a  noise-only  sequence.  This 
occurs  because,  as  can  be  seen  from  Eqns.  (2-17),  a  higher 
value  of  m  lowers  the  upper  threshold  initially,  although  it 
is  then  increased  more  quickly  as  more  samples  are  taken.  A 
preferred  value  has  not  been  decided  upon,   but  a  signal 
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level  equal  to  the  rms  noise  level  or  one-half  this  value 
was  used  for  most  tests. 

Sampling  Rate.  The  sampling  interval  used  in  the  tests 
that  are  presented  was  the  shortest  possible,  being  the 
sample  interval  of  the  discrete  signal  representation,  or 
one  sample  per  cm"  .  Slower  sampling  greatly  reduces  the 
ability  to  detect  low  peaks,  so  the  faster  sampling  rate  was 
maintained. 

Theoretical  Error  Probabilities.  The  intended  probabi- 
lities of  Type  I  and  Type  II  errors,  denoted  a'  and  b' 
respectively,  are  used  in  setting  the  thresholds  as  shown  in 
Eqns.  (2-12)  and  (2-17).  These  equations  show  that  the 
effect  of  a  smaller  a'  (probability  of  false  detection)  is 
to  raise  the  upper  threshold,  without  significantly 
affecting  the  lower  threshold,  thus  making  decisions  in 
favor  of  H,  more  difficult.  The  main  effect  of  a  smaller  b' 
(probability  of  false  dismissal)  is  to  decrease  the  lower 
threshold,  making  dismissals  of  H^  more  difficult.  In  the 
actual  results  the  effect  of  decreasing  a'  seemed  to  be  to 
make  false  detection  less  likely  and  also  to  make  detection 
of  signals  more  difficult.  The  effect  of  changing  b'  was 
not  as  clear. 

4.2  Results  of  Tests 

Complete  sets  of  plots  for  three  different  signal 
sequences,   showing  signal  alone,   signal  with  noise,  filter 
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output,  and  detector  output,  are  presented  in  Figs.  4.1 
through  through  4.22.  The  numbers  printed  along  the  upper 
edges  of  the  plots  are  the  SNRs  for  the  respective  peaks. 
For  each  sequence  several  detector  outputs  are  shown  to 
demonstrate  the  effect  of  changing  some  of  the  detector 
parameters.  In  all  of  the  sequences  the  pulses  are  of  the 
Lorentzian  shape.   The  sequences  are  as  follows: 

(1)  A  set  of  uniform  pulses  with  an  amplitude  of  0.7. 

(2)  A  set  of  randomly  generated  pulses  with  amplitudes 
between  1 . 5  and  3.0.  The  seed  for  random  number 
generation  was  11. 

( 3 )  A  set  of  randomly  generated  pulses  with  amplitudes 
between  1.0  and  2.0.  The  seed  was  11,  following 
the  generation  of  Sequence  ( 2 ) . 

The  amplitudes  referred  to  here  are  set  with  respect  to 
the  rms  noise  level. 

Sequence  ( 1 ) .  This  sequence  was  used  to  test  the 
detector  performance  at  low  SNRs.  The  peaks  vary  in  SNR 
from  0.13  to  0.25.  What  appears  to  be  the  first  peak  in  the 
filter  output  (at  about  80  or  90  cm  )  is  actually  a  noise 
pulse,  as  can  be  seen  from  the  fact  that  a  similar  pulse 
appears  at  the  same  position  in  the  filter  outputs  for  the 
other  two  sequences.  The  actual  signal  peak  follows  it  and 
is  too  low  to  be  detected. 
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Here  peaks  as  low  as  0.2  can  be  detected.   The  second 
and  sixth  peaks,  with  SNRs  of  0.13  and  0.18  respectively, 
cannot  be  detected  with  values  of  a'  anywhere  from  10 
(Figs.   4.5  and  4.6)  to  10   ,   but  they  can  be  detected  when 

a'  =  lO"-*-  (Figs.  4.7  and  4.8). 

-3       -1 
Changing  b'   from  10   to  10     (Figs.   4.5   and  4.6) 

caused  the  secondary  peak  on  the  seventh  pulse  to  be 

detected  separately  from  the  main  peak,   instead  of  being 

missed.   Otherwise  the  choice  of  b'  did  not  appear  to  have  a 

significant  effect. 

When  the  assumed  amplitude  m  was  changed  to  0 . 5  instead 
of  1.0,  with  a'  at  lO'  (Fig.  4.8)  the  sixth  peak  seemed  to 
be  detected  easier,  but  there  was  not  much  effect  on  the 
other  peaks. 

Sequence  ( 2 ) .  Using  a  larger  value  of  a'  generally 
helps  the  detector  performance,  but  in  some  earlier  tests  it 
appeared  that  this  also  made  it  more  likely  that  some  larger 
pulses  would  be  merged  together  in  the  detector  output.  To 
test  for  this  problem  in  the  final  version  of  the  detector, 
a  sequence  with  amplitudes  of  between  1.5  and  3.0  and  SNRs 
of  0.6  to  2.7  was  tested. 

Changing  a'  from  10~  to  10~  did  have  the  effect  of 
causing  less  separation  between  some  of  the  output  pulses, 

but  no  problems  resulted  (Figs.   4.13  and  4.15).    At  a'   = 

-5  -3 

10   ,   a  choice  of  b'  =  10   caused  the  noise  peak  preceding 
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the  first  signal  peak  to  be  detected  separately  (Fig.  4.13), 
while  with  values  of  both  b'  =  lO'  (Fig.  4.14)  and  b'  = 
10  (not  shown),  it  was  merged  with  the  first  peak;  so 
again  the  effect  of  b'  was  hard  to  determine.  Changing  m  to 
0 . 5  had  little  effect  except  that  the  first  peak  was  again 
merged  with  the  noise  peak  (Fig.  4.16). 

The  third  peak,  with  a  SNR  of  0.91,  was  not  detected  in 
any  of  the  tests,  since  it  was  lowered  too  much  by  the  trend 
removal . 

Sequence  ( 3 ) .  This  sequence  was  used  to  test  the 
detector  performance  on  a  varied  sequence  such  as  Sequence 
(2)  but  at  lower  SNRs.  Here  the  SNRs  range  from  0.37  to 
1.50.  The  plots  shown  are  for  a'  =  10~  and  b'  =  lO""^  at 
both  m  =  1.0  and  m  =  0.5  (Figs.  4.20  and  4.21),  as  well  as 
one  for  a'  =  10"^  and  m  =  1.0  (Fig.  4.22).  All  of  the 
pulses  were  detected  when  using  both  values  of  a' ,  although 
value  of  10  results  in  weaker  detection.  There  was  also  a 
false  detection  at  the  beginning  due  to  the  effect  of  the 
trend  removal.  However,  the  noise  peak  which  was  detected 
as  a  signal  in  the  other  two  sequences  was  not  detected 
here. 

Noise-Only  Sequence.  It  was  seen  in  tests  on  a  noise- 
only  sequence  that  false  detections  were  likely  to  result 
from  the  noise  peaks  that  appear  when  the  noise  is  filtered 
(Figs.  4.23  to  4.25).   These  false  detections  appear  even  at 
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very  small  values  of  a' ,  such  as  10  .  Since  the  value  of 
a'  does  not  seem  to  have  much  effect  on  these  false 
detections,  it  may  be  best  to  keep  a'  around  10  or  higher 
to  help  in  the  detection  of  low-level  signals.  In  the  other 
tests,  when  signals  were  present,  these  noise  peaks  usually 
did  not  have  a  significant  effect.  These  results  suggest 
that  there  may  be  a  greater  problem  with  false  detections 
when  peaks  are  farther  apart. 
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TABLE  4.1 

Data  for  Sequence  ( 1 ) 

Pulse 

Center 

Amplitude 

Half -Width  at 

SNR 

Number 

Position 

Half -Height 

1 

51 

0.7 

5.0 

0.146 

2 

151 

0.7 

6.7 

0.127 

3 

251 

0.7 

8.3 

0.240 

4 

351 

0.7 

10.0 

0.236 

5 

451 

0.7 

11.7 

0.188 

6 

551 

0.7 

13.3 

0.180 

7 

651 

0.7 

15.0 

0.246 

8 

751 

0.7 

16.7 

0.167 

9 

851 

0.7 

18.3 

0.182 

10 

951 

0.7 

20.0 

0.233 
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TABLE  4.2 

Data  for  Sequence  ( 2 ) 

Pulse 
Number 

Center 
Position 

Amplitude 

Half-Width  at 
Half-Height 

SNR 

1 

50 

1.45 

5.0 

0.555 

2 

101 

1.63 

5.6 

1.461 

3 

191 

1.45 

9.0 

0.905 

4 

303 

1.94 

14.7 

1.725 

5 

377 

1.57 

10.6 

1.267 

6 

442 

2.15 

7.2 

1.680 

7 

551 

2.55 

7.7 

1.921 

8 

695 

1.49 

16.1 

0.862 

9 

800 

2.77 

19.6 

2.542 

10 

925 

2.53 

8.9 

2.605 

11 

995 

2.58 

11.5 

2.652 
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TABLE  4.3 

Data  for  Sequence  ( 3 ) 

Pulse 
Number 

Center 
Position 

Amplitude 

Half-Width  at 
Half-Height 

SNR 

1 

90 

1.61 

10.6 

1.180 

2 

206 

1.81 

19.7 

.  1.503 

3 

300 

1.56 

16.5 

1.189 

4 

433 

1.72 

10.5 

1.769 

5 

524 

1.52 

10.8 

0.857 

6 

673 

1.45 

15.8 

0.822 

7 

775 

1.01 

14.5 

0.405 

8 

883 

0.98 

5.5 

0.336 

9 

1015 

1.22 

16.8 

0.599 
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V.   CONCLUSION 
5.1  Summary  of  Results 

A  study  has  been  made  of  a  system  for  the  detection  of 
signals  typical  of  those  encountered  in  Raman  spectrocopy. 
The  system  consists  of  a  matched  filter  followed  by  a 
detector  based  on  Wald  testing.  The  detector  has  been 
demonstrated  to  be  capable  of  detecting  peaks  with  SNRs  as 
low  as  0.25  to  0.5.  However,  the  reliability  of  the 
detector  appears  to  be  dependent  on  the  particular  form  of 
the  sequence,  which  is  to  say  that  when  there  are  many  peaks 
close  together  having  varying  amplitudes,  some  may  be  lost 
due  to  the  trend  removal  action,  even  some  of  those  with 
SNRs  as  high  as  1.0  or  2.0. 

It  was  seen  from  the  first  sequence  shown  that  peaks  as 
low  as  0.13  can  be  detected  by  allowing  a  large  value  for 
the  false-detection  probability,  such  as  10~  .  However, 
selecting  large  values  of  a'  will  probably  increase  the 
likelihood  of  false  detections.  Generally  it  appeared  that 
larger  values  of  a'  increase  both  the  probability  of 
detection  and  the  probability  of  false  detections,  as  would 
be  expected,  while  b'  and  the  assumed  signal  level  m  did  not 
have  as  significant  an  effect,  although  a  lower  value  of  m 
did  appear  to  help  detection  in  some  cases. 

In  general,  this  system  may  be  a  useful  addition  to  the 
methods  of  peak  enhancement  that  have  already  been  employed 
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in  Raman  spectroscopy,  in  that  it  provides  for  automatic 
detection  of  spectral  peaks  and  gives  fairly  reliable 
results  at  signal-to-noise  ratios  as  low  as  1.0  or  less. 

5.2  Recommendations  for  Further  Work 

Some  further  testing  of  this  system  would  be  needed  to 
fully  understand  its  capabilities.  Tests  might  be  made  on 
sequences  more  closely  resembling  actual  spectra.  Also,  the 
effect  on  the  performance  when  the  noise  variance  is  greater 
than  the  estimate  should  be  tested. 

The  performance  of  a  detector  based  on  a  multiple- 
sample  Neyman-Pearson  test  might  be  tested  for  comparison. 
Such  a  test  would  be  essentially  the  same  as  the  test  used 
in  the  Wald  detector  to  make  decisions  once  the  maximum 
sample  size  was  reached. 

It  might  be  useful  to  study  the  class  of  detectors 
known  as  nonparametric  or  distribution-free  detectors,  which 
have  been  applied  to  problems  in  radar  and  communications. 
The  characteristic  of  such  detectors  is  that  their 
performance  is  not  highly  dependent  on  the  distribution 
functions  of  the  noise  or  signal-plus-noise.  Such  a 
detector  would  be  useful  when  the  noise  or  signal  statistics 
are  not  well  known.  However,  if  these  statistics  can  be 
estimated  with  some  accuracy,  then  it  would  be  advantageous 
to  use  a  detector,  such  as  that  presented  here,  which  makes 
use  of  this  knowledge.  Reference  [3]  is  one  source  for 
information  on  nonparametric  detection. 
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riFTFCTnR..TF:ST.F(1R 
KANBftB    STATE    IJN.T VFRBITY 


C 

C  PROGRAM    riETECTOR_TEST 

C 

c 
c 
c 
c 
c 
c 
c 
c 

c 

PURPOSE!  M.T.in  routine  to  simp  r,-?(.o  risni^l  dnt3»  c^rr'.i  out  filtrrind 
nni:1  the;  detection  r-rorrjduro;  writr.'  d,-jta  to  r,nd  ro.3d  dntn  from  diski 
and  plot  di^ta  ^ccuonrcr. . 


VAX-11  FORTRAN  SOURCE  FILENAMK:: 

DEPARTMENT  OF  EFCF 

author:   n.imon  Mick 

DATE  created:   July  19R7  (Fin.il  vcrr.ion) 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

program  detector.test 
implicit  none 


ROUTINE(S)  CALLED  OR  ACCESSED  BY  THIS  ROUTINE: 

DETECT 

FILTER 

PLOT 

READ-DATA 

RMS..CALC 

SIGNAL_GEN 

URITE_DATA 


REAL      SIGNAL(10?1).  N0ISE(1024)»  PTGN0I(1034)r  OLITPUT  ( 10;.M  )  , 
Y_DATA(  1021)  7  MEAN<1A)7  AMPI..(t6)f  HUHH(IA). 
SNR(16)r  RMS  I  RMS_2»  AMAXr  ALPHA  r  KiETAr  LEVEL 


INTEGER   DET(1021)j  OPER;  SE0..TYPi;:7  SHAPE;  SEED7  ARRAY_NUM. 
I>  Tr  DFVICE_NUM 

CHARACTER*!  A 
CHARACTER*20  FJLE..NAME 

Main  men  u . 

DO  WHILE  (0PFR.NE.6) 

TYPE  «j  'ENTER  NUMBER:  ' 

TYPE  *>  '  <1)  Gcncrste  3    !:;i£!n3lr  .-jdd  noi^ct  end  filter' 

TYPE  *r  '     (?)    Rc?3d  .-5  filterod  r>iannl  <3SQi.rcncc  from  3  file' 

TYPE  *r  '     (3)  Set  threshold  and  rerfonn  detection' 

TYPE  *7  '  (1)  Plot  d.it.^' 

TYPE  *»  '  (5)  Write  n  d.-jta  ccnuenee  to  ,1  file' 

TYPE  *7  '  <A)  Exit  the  rrortrsm' 

READ  *f  OPER 


IF  (OPER.EQ.l)  THEN 

Resd  noise  cenucnce  i?nd  ohtsin  me.-jr>Lirement  of  rmr>  noi!;c  level, 
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c 

c- 

c 


c 

c- 

c 

c 

c- 

c 


TYF'E*»     'Rcsdintf    noise    file.' 
Cf\LL    RFr^ti.n.'^lYK NOISE) 
CALL    RMS_CftLC(K'niSFr     RMS) 

Get    .i  nf  o  rin.Tt  i  on    from    the    u!ser    for    ih(?    dcncr.^tion    of    rulrjcr. . 

TYPE*.     'Enter    <:;h,-5r'G'    of    PLi:tc(r<>;     (1)     Raudcicn       (2)     Lorcntzinn ' 

READ«j     SHi^PE 

TYPE*»     'Enter!     (!)     R.mdom    r'iilr>(vr>       <  2 )     Uniform    ru.ln(\-5' 

READ* 7     SEQ-TYPE 

TYPE*r     'Enter    m  3  >;  i  m  u  m    cimpl  i  tucio    for    rultses!' 

REAn«7     i^HAX 

IF     (SEQ-TYPE, EQ.l)     THEN 

T  Y  F'  E  *  >     'Enter    r.  p  G  d    for    ?■  i  c!  n ."!  1     f.(  p  n  e  r  o  t  i  o  n  !  ' 

REAIi*,     SEED 
END    IF 

CALL    SIGNAL..GEN(SIGNALf     NOISFr     BIIAPF;     SFn.TYPEr     AMAX» 
RMS,     MEAN,     AMPL,     HUHH>     SNR7     SEED) 

no       1^1  r     1024 

SIGNOKI)     =    SIPNAL(I)     +    NDISEd) 

ENC  no 

Obtain  noicc  vsrisncc  r.t  filt(-:r  outr-i.it  for  ur.p  in  detection  routine 

CALL  FILTER  (NOISE.  OUTPUT) 
CALL  RMS..CALC(OUTPUT;  RMS_2) 

Obtain  o u  t r- u t  of  matched  f  i  1 1 c  i- . 

CALL  FILTER(SI6N0Ir  OUTPUT) 

Write  r-iilGG  data  to  r.    file. 

OPEN  (UNIT-lf  STATUS-'MIIW  r  T  ILE  •  '  PUI  SFS  .  DAT  '  ) 
DO  I  ^  1;  15 

IF  (MEftN(I)  .NF.O.  )  THEN 

URITEdi     1)     Ir     MEAN(l)t     AMPL(I)t     HUHH(T)i     BMR(I) 
FORMAT     (11,27     IFIO.S) 
ENIi    IF 

END  no 

WRITE     (Ir     1)     lAf     0,0>     0.0.     0,0t     0,0 


C 

c- 
c- 
c 


ENn     IF 

IF     (0PER.EQ,2)     THEN 

Noise    seoencc    tg    needed    to    c."<lcu];:te    estimate    of    noiriC 
standard    deviation. 

TYPE*»     'Recuire    noir^p    seouence    ,:>nd    r>irJria]     t>/pc    u!5Pd    in    ' 

'orisinal    <;ianal  .  ' 
CALL    REAn_nATA(NOISE) 
CALL    FILTER(N0ISE7    OUTPUT) 
CALL    RMS_CALC( OUTPUT T    RMS. 2) 

Obtain    file    of    pulse    data    r>ri    thai.    SHHr,    can    be    plotted, 

TYPE*7     'Enter    filename    for    rul<ie    dat,?.' 
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c 

c~ 

c 


READO,  3)  FILE_Hi^ME 
FORMAT  (A20) 

0PEN<UNIT=1.  STATUS^' 01.11' I  F  II.E-FII.E  ..NAMF  ) 

no  WHILE  (I  .NE.IA) 

REAIKl,!)  I,     MFAN(I)i  AMPl.  (  I  )  ,  IIWHH(I)>  SNR(I) 

ENn  no 

Read  file  of  filtprcd  r.iSnsl. 

CALL  REAn_nATA( OUTPUT) 
END  IF 

IF  (OPER.EQ.?)  THFN 

TYPE*?  'Enter  desired  error  r  rohrsbi  1  i  tics  (Type  Ii  Thft  II) 

REAn*;  ALPHA;  BETA 

TYPE*»  'Enter  assumed  si^nol  Ir^vpl!' 

READ«»  LEVEL 

TYPE*»  'Enter  s.-jmr-l  ins  interv.il!' 

REAn*>  T 

CALL  nETECT(OUTPUTr  DFT,  RMS_2f  ALPHA.  BFTAi  LFyFL r  T) 
ENn  IF 


IF  ( (0PER.EQ.4) .OR. (0PER.F0.5) >  THEN 

15         TYPE*,  'Select  dntn  3rrn\it' 

T  Y  P  F  *  I  '  (  1  )  K  i  d  ri  a  1  without  n  o  i  r.  c  ' 

TYPE*7  '  <;.')  Signal  with  noi-5f?' 

TYPE*.  '  (3)  Filter  outrut' 

TYPE*.  '  <<1)  nritector  output' 

REAn*.  ARRAY..NUM 
ENi:i  IF 

IF  (OPER.EQ. '!)  THEN 

IF  (ARRAY-MUM.EO.l )  THEN 
no  I  =  1 .  1024 

Y_nATA(I)-^SIGNAL(I) 

ENn  no 

ELBE  IF  (ARRAY_NUM.F(3.2)  THEN 
no  I  -  1.  1034 

Y_nATA(I)^SIGNni<I) 
ENn  DO 

ELSE  IF  (ARRAY_NUM.En.3)  THEN 
no  I  -  1.  1024 

Y-nATA(I)=OUTPUT(I) 

END  no 

ELSE  IF  <ARRAY_NUM.EQ.4)  THEN 
no  I  -  1,  1024 

Y_riATA(I)-nET<I) 

ENn  no 

ENn  IF 

TYPE*.     'Fntcr    dpvirc    nunihert     (1)     7'!75       (2)     'lOll' 

READ*.     nEVICE-NUM 

CALL    PLOT(Y_nATA,     nryiCE_NUM.     AKRAY..NUM,     MEAN.     SNR) 
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TYPE  *7  'Ho  'jou  want  to  mnke  .'jnothsr  riot,?' 
READ(5>  2)  A 
FORMAT  (.11) 

IF  (A.EQ.  'Y'  .■>  RO  TO  15 
END  IF 


IF  (OPER.EQ.fi)  THEN 

IF  <ARRAY_NUM.E0.1)  THEN 
CALL  URITF_.DftTft(S:((5N(^,l.  ) 

(ARRAY..NUM.E0.2)  THEN 
URITE_r.iATA(SIGNni  ) 
(  ARRAY.. NUM.  FQ.  3)  THEN 
lJRITE_nATA(  OUTPUT) 
(ARRAY...NUM.EQ.-1)  THEN 
WRITF-IiftTAdiFT) 


ENIi 


ELSE  IF 
CALL 

ELSE  IF 
CALL 

ELSE  IF 
CALL 

ENIl  IF 

IF 


END  DO 

STOP 
END 
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c*«***«**«*««*«»******««**«««««**r.^^f*:********r***««**«***«*«**«*^ 


SUBROUTINE  REAti-nftTA 

yAX-.tl  FORTRAN  SOURCE  FII.  ENAMh'. 

DEPARTMENT  OF  F.ECF 

AUTHOR!   ricmon  Mick 

DATE  CREATED!  .hi\'.;    1987  (Firml  verr.ion) 


RlJ._nATA.FnR 

KANSAS  STATE  UNIVERSITY 


PURPOSE!   This  routine  rcndr,  .■;  drtn    r,(i  rurncc  from  r,    di!:^K  file. 

ROUTINE(S)  CALLED  OR  ACCESSED  BY  THIS  ROUTINE! 
DISK  IN 


real 


C****«******«***«*********«««*«****«««*******«****«****************^ 

C 

C  CALLING    SEOUENCE!        CALL.    READ-DATA     <SEI.niEMCE) 

C 

c 
c 
c 
c 
c 
c 
c 
c 

SUBROUTINE  RE AD_ DAT A ( SEOUENCE ) 
REAL  SEQUENCEdOSI) 

CALL  DJSKIN( 'REAL' r  SEOUENCE,  10?1j  LHNACT,  'PROMPT' r  'PROMPT', 
g   I ERR  7  0) 

RETURN 
END 


ARGUMENTS  REQUIRED! 

SEQUENCE         output 
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RU_riATrt.FnR 

KAMSriS    STATE    DNiyFRSTTY 


C 

C  subroutine:    URITE-TiATA 

c 
c 
c 
c 
c 
c 
c 
c 
c 

c 

CALLING  SEQUENCE!   CALL  URITE..riATA  ( SLRUFNCE ) 


yAX-11  FORTRAN  SOURCE  FILENAME! 

DEPARTMENT  OF  EECE 

AUTHOR!   riamon  Mick 

DATE  CREATED!   .h.il^i     1987  (Fin,-;!  vpirion) 


PURPOSE!   This  routine  uritcc  .:;  cJrt-,-!  TTCciicnce  to  rs  disk  file. 

ROUTINE(S)  CALLED  OR  ACCESSED  BY  THIS  RnilTINF! 
DISKOUT 


input      real 

A  r  r  >T.i  f  r  CI  Hi  which  data  i  r  to  bo  written. 


ARGUMENTS  REQUIRED: 

SEQUENCE 

C#****.***********»«***.*«****  **«*«*  *****«*«***********«««««««***« 

SUBROUTINE  WR I TE  DATA < SEQUENCE  ) 

REAL  SEOUENCE(1021) 

CALL  niSKOU( 'REAL' !  SEQUENCE.  1024r  'PROMPT' r  'PROMPT', 
S.       I  ERR,  0) 

RETURN 
END 


62 


c 

SUBROUTINE  SIRNAL-RfTN 


VAX- 11  FORTRAN  BOURCF  FILFNAMFT! 

DEPARTMENT  OF  EFCF 

AUTHOR!   Damon  Mick 

DATE  created:   July  1987  (Final  vcir.ion) 


BJGNAI._(5FN.FnR 

KANSAS  STATE  UNIVERSITY 


C 

CALLING  SEQUENCE!   CALL  SIGNAI-CFN  (SIGr  NOIt  SHftPEr  SFG_TYPFr  A«AX» 
RMS  7  MEAN  J  AMPLj  HUHH;  SNRi  IX) 


PURPOSE!   Thic  routine  drncrotcs  r.prici;  of  rultscr.  hrvinii    cither 
Lorpntzinn  or  G.iur>r>i.Tn  ■sharti'j. 

ROUTINE(S)  CALLED  OR  ACCESSED  BY  THIS  ROUTINE! 
RAMDU 


ARGUMENTS  REQUIRED! 
SI6 


NOI 


SHAPE 


SEQ-TYPE 


AMAX 


RMS 


MEAN 


AMPL 


HUHH 


SNR 


output,     real 

Tho  nrr.T.f  contsinina  tho  r.ian.Tl  r>Gauoncer.!  without 

noiso  iiiddcd. 

i  n  r  u  t      real 

T  h !?  n  r  r  .T.)  c  o  n  t  n  i  n  i  n  "3  t  h  o  n  o  i  o  o  s  o  a  u  c  n  r  e  . 

input      i  n  t  G  ?5  (T  r 

Indirnt:?r>    the    r>hnr'c?    of    the    rul-ser.    to    hn    -.ir-no  rat!?d . 

1!       R.iucriirn    pulr>pc. 

2!       Lorontzi.Tn    Fulcss. 

input  intodpr 

Indicator)    t h :?    t 'j r  c!    of    T)  t?  r  i  i? ■-.    to    bo  3 c n o  r  a  t o d  . 

1!       Rnndortili'    dcnoiatird    pulr.cs. 

'<!  I       Pulr>f;r>    uniform    in    hoilht    snd  spscins. 

input  rejil 

The    ra.T>!imum    nmplitudj;    for    the    pul'^ss. 

i  n  p  LI  t  real 

The    rnir>    noi^.^    level. 

output  rpsl 

Th:?    nrr.T.i    of    tho    ccntc?r    por. itions    of    nil    tho    pu15e?5> 

o  u  t  p  u  t  r  <?  ii  I 

Tho  nrrni!  of  tho  nmp  1  i  tudor.  of  all  the  Fulr.c;. 

output     real 

Tho  <Trr3'.(  of  tha  huhh'r.  of  all  tho  pulsos. 

output     r  p  5 1 

Tho    arra>.f    of    tho    SNRr.    of    all    tho    pulr.os. 
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C  IX  input      intcSer 

C  The?  r.t^.T'd  far  r.mdom  numbrr  aono  ration. 

C 

SUBROUTINE    SIGM.'il.  ..IjENC  St  fi  j     NOI;     SHAPIf:;     SEQ.TYPI?: ;     .'iMAXt     RMS. 
-       MEAN.     AMPL.     HUHHt     SMR.     .TX) 

IMPLICIT  NONE 

REAL      SIG(10?'»).  Nni(102'1).  BGPEAK(IS),  MFftN<«). 
AMPL<»)7  HUHH(«)»  SMR(»).  AMAXi  YFLj  RMS. 
SSUM.  NSUMi  SHEW 

INTEGER   SEO-TYPE.  SHAPE;  IX. 

lY,  PNUM.  I.  J.  MIN-POINTr  MAX...P()INT 
PRINT*.  IX 
C 

C —    Removp  nnu  p  rcviour.l'j  rttrrie  i-istcd  !ii3n.Tl  from  t.hp  nrrs-j 
C--    SIR  oncl  3ns<  pulr>c?  dnto  from  tha  riilr.f?  datn  am>\i<i, 
C 

DO  I  ^  1.  iO?.A 
SIG<I)  ^  0.0 

END  no 

PNUM  =  0 

DO   I  -  1.  15 

SGPEAK(I)^0. 

MEAN< I)-0. 

HUHH( I)=0. 

AMPL<I)^0. 
END  DO 

IF  (?.FO_TYPE.EC).l  )  THEN 
C 

C      Procedure  for  r ."!  n  d  o  id  r-  u  .1  T)  c  s  . 
C 

C         Obtain  position  of  fir;;t  puItjC?. 
C 

CALL  RANDU<  IX.  lY.  YFI  ) 

MEAN(l)  ^  50  +  100*YFL 
C 

C  Obtain    half    width    ijt    half    hci.dht    for    current    rulce    and 

C  calculatf?    ^t ,    dcv .     from    hwhh. 

C  Obtain    Pulr>G    hciaht.     uith     rer,r(?ct    to     riiir>    tioir.G     lovcrl. 

C 

DO    UHILE     (PNUM, LT. 15     .AND.     MEAN  (  PMUM+1  ).  LT  .  .1  024  ) 

PNUM    -'    PNUM    +     1 

CALL  RANDIKIX.  lY.  YFI.) 

HUHH<PNUM)  -  5.0  +  15.0*YFL 

CALL  RANDUdX.  lY.  YFL) 

AMPL  (PNUM)  -■     (  1  .+YFL)*0.5*AMAX»RMS 
C 

C  Generate  the  point  values. 

C 

MIN..POINT  -  MAXd.  N  INT  (  MEAN  (  PHUM  )  )  -  1  00  ) 

MAX...POINT  ^  MINd.O?'!.  NIMT  (  MEAN  (  PNUM  ))  f  100  ) 


en 


no  J  -^  MIN._Pnt«Ti  MAX.POINT 

IF  (SHAPE. FQ.t)  THEN 

G  3  u  s  G  i  n  n  r  u  ]  5  c  ?  • 

SDEV  ^  HUHH(PNIJM)/SPRT(2.*L0R(?.  )  ) 

SIG(J)  -^  SXR<,))+AMPL.(PNUM)«EXP(-<  J-MEAM(PNUM)  )««2/ 
(2*StiEy*«2)  ) 


ELSE  IF  <BHAPF.ER.2)  THEN 

Lorentninn  pulr>Gr>. 

SIG(J)     ^    SIB(J)+AMPL(PhMm)/ 

(l  +  (  (  J-MEAH<PMl)M)  )/HUHH<PNIJM)  )««?) 

END    IF 

END  no 

Obtain  the  rar.ition  of  thp  n(?;;t  fiiilrc. 

CALL  RANnil  (IXt  IY>  YFL) 

MEAN(PNUM  +  1)  ^  MEAMCPNUM)  +  50  +  10()«YFL 

END  no 

ELSE  IF  (SE0..TYPE.EC1.2)  THEN 

Generstion  of  uniforin  r-ulncG. 

no  PNUM  :^  Ir  10 

MEAN(PNUh)     ^    51    +     (PNUM-.t  )«100 
HUHH(PMUM)     ^    5.0     ^     (  PMUM--1  )  *15  .  0/9  .  0 
AMPI.  (PNIIM)     ^    AMAX 
SnEV    -=    HUHH<PMUM)/SQRT(2.«L0G<2.  )  ) 


DO    J 


(PNUM-l)«100  +  lr     PK'tlH*100 


IF     (  SHAPE.  ECl.l)     I  MEN 

SIG(  J)     '■■    AMAX«FXP(-(.I-MFAM(PM1IH)  )««?./(  2  .  *f>ni:y««2  .  )  ) 
ELSE 

SIG<J)     -    AMAX/d  .  +  (  <  J-MEAN(PNIJM)  )/HlJHH(PNIJM)  )««2.  ) 
END    IF 

END    DO 

END  no 

PNUM    =10 
END    IF 

IF     (SHAPE. EQ.l)     THEN 

TYPE    *»     'BiariHl    dcnerjition    (G.-jusGi.in)    cnmrlptGd.' 
ELSE    IF    (SHAPE. EQ. 2)     THEN 

TYPE    *  r     '  S  i  r)  n  3 1    S  p  n  c  r  n  t  i  o  n    ( I.  o  r  c  n  t-  z  i  .i  n  )    c  o  in  r-  ]  r  t  p  d  ' 
END    IF 
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Calculate  SNR  for  bII  r-ulspr.. 

no  I  -  1,  PNUM 
SSUM  -^  0 
NSUM  ^    0 

CO  J  ^    NINT(Mh:AN(I)-2«HUHH<I)  )  ,  MINT  (  Mf  AN  (  I  ) +2*HUHH  (  1  )  ) 
IF  (  (J.CiE.n  .AMn.  <  J.L.E.10;3'1)  )  THFN 


SSUM  -  SIR(J)«*2 
NSUM  -  N0.T(J)«*2 
END  IF 
ENn  DO 

SNR(I)  ^    SSUM/NSUM 
ENi;r  DO 
RETURN 
END 


SSUM 

M!:>IJM 
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c 

C  SUBROllTINfT    RftNDU 

c 

C        yAX-l.t  FORTRAN  StlURCF  FIIENAMK!  SIRNAI.  .GF  N  .  FOR 

C        DEPARTMENT  OF  EECE  KANSAS  STATE  UNiyFRSITY 

C 

c 

c*««*«***«*««**«»«*«««««««***«*«*««***«««***««»***«****«**************** 
c 

C        CALLING  sequence:   CALL  RANnU  (IX>  lYr  YFI.  ) 

c 

C        purpose;   ThiG  routine  acncratec  .i  number  uhich  h.-?r,  n    v.tXuc    evenly 
C        dir>t,rihut,E?d  betuopn  0  ond  1. 

C 

C        ROUTINE(S)  CALLED  OR  ACCESSED  EY  THIS  ROUTINE: 

C  None. 

C 

C        ARGUMENTS  REQUIRED: 

C 

C  IX  inrut  integer 

C  Tho    r.ord    for    the    random    ni.imbor    l^nerstion. 

C 

C  lY  output  inteaer 

C  Th!7    rnndom    numbrr    thnt    i '-.    flcnoratt^d?    boforo    b:?in3 

C  normclized. 

C 

C  YFL  outrut  integer 

C  lY    normnlizsd    to    s    v.tIm:?    bf^twesn    0    and    1. 

C 

C**#«««*«*««*«*«*««*«««»«««««*******««««*««««***«« ««««*««*«««««« 

SUBROUTINE  RANDU  (IX 7  IY7  YFL) 

C 

C —    Random  number  flenerntion  routine.   Vnlue  of  YFL  will  be  evenly 

C —    distributed  between  0  nnd  1. 

C 

lY  -  IX«65539 

IF  (lY)  5f  6  J  6 

5  lY  =^  lY  +  21')7'183A^7  +  1 

6  YFL  ^  lYmO.IASdAlSE-? 
RETURN 

END 
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WAX- 11  FORTRAN  SOURCF  FILFNAMF: 

DEPARTMENT  OF  FECF 

author:   Iismon  Mick 

DATE  created:   Ji.ila  1987  (Final  vpir.imi) 


FILTER. FOR 

KANSAS    STl^TE    1JN.TVFRSITY 


C 

C        SUBROUTINE  FILTER 

C 

c 
c 
c 
c 
c 
c 
c 
c 

c 

C        CALLING  sequence:   CALL  F IL TER < SIO  i  OUTF) 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c*****« ****«*«**«««***«««»««*«#* ********** ***««*«*«*«««««**««««* 
SUBROUTINE  FILTER  <SIG»  OUTP) 


purpose:   Thi?  routine  computes  the  rri-riuli-  of  f'i  1  tc  ri  nrJ  the  di^ta 
saaucnco  SIG  uith  n  filter  m.3tchc;r|  to  s    Lorpntzian  pulse  ch.TPe. 

ROUTINE(S)  CALLED  OR  ACCESSFn  BY  TKIB  ROUTINE: 
FFT 


ARGUMENTS  REQUIRED: 
SIG 


OUTP 


input      real 

The  nrr.T.f  of  the  data  to  be  filtered. 


output     reol 

The  srr.TM  contninins  the  output  of  the  filter. 


REAL      SIG(1021).  0UTP(1021)<  FIL(i;>8)i  niH-y 
COMPLEX   H( 256)7  X(256)7  Y(2SA) 
INTEGER   J.  INU 


G e n e r fs t  i o n  of  filter  r e r. p o n •:- c  . 

DO   J  =  1.  128 

FIL(J)  =  1/(1  +  ( ( J-63)/10.0)*«2) 
END  DO 

Initialize  the  filter  output  orra-f. 

DO  J  =  If  1021 
OUTP(J)  ==  0 
END  DO 

Transfer  filter  rcGPonnc  to  coniple;;  H  orriT^i  and  ini  ti  ril  i;::r 
second  h-ilf  of  H  array. 

DO   J  =^  Ir  128 

H(J-)  rr    CMPLX(FIL(  J)  .  0.0) 

H(.J  +  128)  ^  (07  0) 
END  DO 

Convolve  filter  •scQucnce  with  sisnal  Eenucnr.e  b'.i    ovcrl  ap-rsdd . 
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c 
c 
c 


Seancnt.r.    of    r>i3ri.il    r>r;aucncp    arc    placed    in    comrlor    X    3rr3u. 

INV=0 

CALL    FFT    (H.    256.     INV) 

no  Q  -  07  102-1/128-1 
Do  the  Qth  Gcament. 

DO   J  -^  1-  128 

K  ^    ,J  +  Q*128 

X(J)  =  CMPLX(SIG(K)  r  0.0) 

X<J+128)  -  (0;0) 
END  DO 

IM«^  =  0 

CALI  FFT(X.  256j  TNV) 

X  and  H  orro'jr,    contnin  the  DFTr  of  the  drjtn  nc.dmcnt  .md  the  filter 
reGrons3»  rospcctivf!l>j, 

DO  J  ^  1 r  256 

Y<J)  ^    X(J)*H<J) 
END  DO 

Take  inverse  tren^form  of  Y. 

INV-1 

CALL  FFT(Y.  25Ar  ItW) 

DO  J  ^  1 7  256 

K  =  J+Q*128 

OUTP(K)  -  nUTF(K)+RFAL( Y< J)  ) 
END  DO 


END  DO 

RETURN 
END 
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c 

C        SUBROUTINE  DETECT 

C 

C        VAX-11  FORTRAN  SOURCE  FILENAME:  HETECT.FOR 

c 

C        DEPARTMENT  OF  EECE  KANSAS  STATE  UNIVERSITY 

C 

C        author:   Damon  Mick 

C 

C        DATE  created:   Julw  1987  (Final  version) 
C 

C*«*«****»*«««*««***«*««#«#******t.*«««*«««««**##**«««*»*****»«««« 

c 

C         CALLING  sequence:   DETECT  (  OUTP  ,  DFTr  RMKi  ALPHA.  JHKTA.  lE'-lFL, 

C        SAMPLE_INTERyAL) 

C 

C  purpose;       This    routine    pprformr    n    Urjilri    trrtina    rrocedurp    on    l.hc    ouir-ut 

C  of    tho    matched    filter. 

C 

C        ROUTINE(S)  CALLED  OR  ACCESSED  PY  THIS  ROUTINE: 

C  None  • 

C 

C        ARGUMENTS  REQUIRED: 

C 

C  OUTP  inr-ut  resl 

C  Tht?    <Trr.Tvi    cont.Tinina    tho    matchGd    filter    outr-ut. 

C 

C  DET  outnit  intcctcr 

C  The    nrr.T.)    cont.iinincf    thj?    detector    output. 

C 

C  RMS  i  n  r  u  t  real 

C  The  rmu  noir.e  level  nt  the  m^tched  filter  output. 

C 

C  ALPHA  input      real 

C  The  theoreticnl  value  for  the  f  al  se-.Tlorm  prohah  i  I  i  t«  . 

C 

C  BETA  input      real 

C  The  theoretical  value  for  the  f a  1  co-di smisssl 

C  p  rohahi  1  i  tL< . 

C 

C  LEVEL  input      real 

C  The  ar.r.i.imed  value  for  the  <5i3nal  level. 

C 

C  SAMPLE_INTERVftl.     input  intescr 

C  The    r.pacin.r!    hetuesn    rsample-;. 

C 

SUBROUTINE  DETECT(OUTPi  DET7  RMS,  ALPHA,  BETA,  LEVELf 
-   SAMPLE-INTERVAL) 

IMPLICIT  NONE 

REAL  0UTP(102'1)t  RMS.  ALPHA,  BETA.  I  FVFL  r  JWIT_BIAS_0, 

8  INTT..BIAS..1  ,  BIAS..O,  BIAS.).,  SUM,  EXPT-VALUE 

INTEGER  DETdOS-l),  SAMPLE,  INTERVAL  r  SAMPLE...NI)  r  DI;TECT_FLG, 

S  I,  J,  K,  MAX..SAMPLE...LEN 

HAX_SAMPLE_LEN  -  .15 
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c 

C      Hetermination  of  thrcGhold  levels. 

C 

INIT-FiIrtS.O    ^    RMS*L0G(BETA/(1  .-Al  FH^;)  ) /I.EyEI  +0  .  5*LFVEL*RMS 
INIT-BIhS...l     ^    RMS*LOr><  (  1  . -BETA  ) /ALPHA  ) /LEyPL  +  O  .  5*LEVEL*RMS 

C 

C  Procedure    for    firct    decision. 

C 

I    ^    1 

DETECT-FLG    =    0 

SAHPl.E-NO    ^    1 

BIAS_0  -^  INIT_BIAS_0 

BIAS_1  =  INtT_BIAS_l 

SUM  =  0. 

no  WHILE  (DETECT_FLG.EQ.O) 

SUM  ^  SUM  +  OUTP(I) 
C 

C         Decisions  made. 
C 

IF  <SUM.LE,BIAS_0)  THEN 
C 

C  No  sisnsl  r-resent. 

C 

BETCI)  -  0 
CiETECT_FLG  -  1 

ELSE  IF  (SUM.GF.BIAS..1)  THEN 
C 

C  S  i  a  n  3 1  r-  r  e  !i  c  n  t , 

C 

BO  K  -  I-SAMPLE.-INTERyAL  +  ]  r  I 

DETCh)  ^  1 
END  DO 
DETECT-FLG  -  1 

ELSE 
C 

C  No    dociGionr^    made!    ui>d,it,G    thrcr,holds    ,7nci    r.rjmr-lc    rount. 

C 

DET(I)    ^    0 

SAMPLE.NO  =  SAMPLE..N0  +  1 

BIA5.0  -  BIAS..0  +  0.5*RMS*LEyEL 

BIAS-l  =  BIAS_1  +  O.S«RMS*LEyf:L 

END  IF 

I  =-  I+SAMPLE_INTERVAL 

END  DO 

no  UHILE  (I.LE.1024) 

DETFCT^FLG    ^    0 
BAMPLE_NO    :^    1 
BIAS..0    ^    TNIT.BIAS_0 
BIAS.l     =    INIT_.BIftB_l 
SUM    ^    0. 
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c 
c 
c 

c 
c 
c 


no  WHILE  (  (DirTR-CT.Ft.G.Pn.O)  .ftND.  (I  .I.E.  to;M) ) 
SUM  -■    SllM  +  OUTP(  t) 

IF  (SUM.LF.BiriB.O)  THEN 

nFTECT_FLG  ^  1 

no  J  =  I-SAMPI.F_INrFRyftl..«SftMPI.E_K"0+lr  I 
DET(J)  -  0 

END  no 

PRINT*,  nFT(I),  SAMPLE-NO 

ELSE  IF  (SIJM.(iE.BIAS-t  )  THEN 

DETECT-FLO  ^  1 

no  K  =  I-SAMPLE-JNTERVAt«SftMPLF_NO  +  l  ,  I 
nET(K)  ^  1 

END  no 

PRINT*.  nFT(T)i  SAMPLE_NO 
ELSE 
No  dpcision  mnde . 

IF  <SAMPLE_NO.LT.MAX..SAMPLE.J.FN)  THHN 

Continue    G.Tmr-l  ind !    urdatc    thrci;hnldr>    ,-!nd    cnmr-lc    count. 

SAMPLE- NO    -^    RAMPLF..Nn  +  l 
BIAS..0    -    BIAS..0    +    0.5«LFyFL*RMS 
BIAS_1     ^    BIAS-I     +    O.S«LFyFL*RMS 
ELSE 

Maximum    £.-?mr  ]  e    Icnsth    rp,-!chcdf     m.ikc    .n    dpcir.jon    h.-3!SGd    on    3 
Noiiraon-Poorson    tr-rit    for    Gamrlo    •r.ir.'^    MAX  ..SAMPLE... LEN . 

IiETECT.FI.  (3    ^    1 

IF  (RUM.GF.  <3.P.«SQRT(FI  OAT<MA>;..SAHFI  B...  LFN)  )«RMS)  )  THEN 
DO  J  --     r--HAMPLF.INTERVAL*SAM|-'LE...NO+l  7  I 
nET(J)  -  1 

END  no 

ELSE 

DO  J  -  I-5AMPI.  F_IN1ERVAI.*SAMP1.E_N0  +  1  r  I 
DET<J)  ^  0 

END  no 

END  IF 

PRINT*.  DFTd) 
END  TF 

END  TF 

L o c £5 1 i  o n  of  n p ; ; t  r> .■j m r  1  o  . 

I  =    I  +  SAMPLE...TNTERVAL 

END  DO 

END  no 

RETURN 
END 
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c****«*****«* **************************************** ******** 

c 

C        SUBROUTINE  PLOT 

C 

C  VAX-11    FORTRAN    BOURCF    F.TLENr'.MFr !  PinT.FOR 

C 

C        DEPARTMENT  OF  FECE  Kf-.KSf-.S    ST(^TF  UWCVFRSITY 

C 

C  AUTHOR!        Di5inon    Mirk 

C 

C  DATE    CREATEfi:        July    1987     (Fini?!    vprr.ion) 

C 

C****«*««*«*««**«******««««««««** ********* **««*««**««*««******************* 

c 

C         CALLING  sequence:   CALL  PI.  (IT  ( Y..I.V'.Ti*;  r  JiFV  .NUMBER  r  PI  OT.  NUMBER. 

C        MEAN  7  SNR) 

C 

C  PURPOSE!       This    routine    rlotc    dn.tB    neGuerirc<7.    rcrr  rcGcciti  n:!    pither 

C  nian.il;    sidn.il    rlur,    noino;    niGtchori    filt.or    outruti    or    dotoctor 

C  outputj    35    chosen    in    the    mcin    r  rorfrrsm . 

C 

C         ROUTINE<S)  CALLED  OR  ArCFSSFIi  BY  THIS  RtUJTINF! 

C 

C  GETUTX 

C  PAX  IS 

C  PCLOSP 

C  PCLRSC 

C  PFIXNM 

C  PFLTNM 

C  P I N I T 

C  PLINE 

C  P0RI6 

C  PPLOT 

C  PPLOTR 

C  PSCALE 

C  PSLPEN 

C 

C        ARGUMENTS  RFQUIREn! 

C 

C  Y_DATA  inrut      real 

C  Arr.-!>.f  contnini  nl  thr>  drstn  to  he  plotted. 

C 

C  DEV-NUMPER       input      integer 

C  Indicntcr.  the;  rlottina  dovice, 

C  1!      717fi  Plotter 

C  2!      401'1  nisrl.iy 

C 

C  PLOT_NUMBER      input      integer 

C  Indionti;;r>  tho  d-itn  uhich  in  bcina  rlottod. 

C  1  1      Sian.Tl  nlone 

C  2!      5ir:!ni-?l  r-lus  noi  =  =' 

C  3!  Matched    filter    outr-ut 

C  4  J  Dotf?ctor    out!  lit 

C 

C  MEAN  input  real 

C  Arr.T.!    cont.iinina    thp    c;7ntf?r    rocitions    for    o.ich 

C  pulse . 

C 

C  SNR  inrut  real 

C  A  r  r  n '.(    c  o  n  t  <i  i  n  i  n  3    t  h  o    s  i  fl  n  o  1  - 1  o  -  n  o  i  s  3    r  -3 1  i  o  -;    for 
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C  G  a  c  h  p  u  1  s  c  . 

C 

C***i|c**#«««»!!t«««)|t*«»**#«*«*«*«*#«  ««**««******«:  «*«*«***#« 

SUBROUTINE  PLOT  (  Y..riATA  j  lUf.y  ..NUMBER  >  FL.OT...NUMBEK  ;  MEAN.  SNR) 

IMPLICIT  NONE 

REAL       X_IlATA(  1021  )  r  Y-D  AT  A  <  1  021  )  t  Mi:.AN(«)r  SNR(*) 
INTEGER   i;iEV...NUMBER7  PHH  ..NUMBER  j  inOT7  I 

REAL      X,  Yr  FIRSTX.  DELTAX.  lUVLNXr  KIRSTYr  DFLTAY, 
niVLNY7  FTRIiEI.<  1)  7  SCALE 

CHARACTER«20   STRj  UNITS 
CHARACTER*f)    INARY 

HO  I  ^  1.  1024 

X-DATAd)  ^  I 
END  DO 

IF  (hev-number.eq.i)  then 

DEy-NUMBER  ^  7-175 
SCALE  -^1.0 

ELSE 

DEy_NtJMBFR    -    1014 
SCALE    -    1,52 

END    IF 
C 

C--   Initialize  rl  ot  tir-i.c!  device'  cet  orin.tn. 
C 

CALL  PINIT(nEV_NUMEiER.  '  't  SCALEr  'ft') 

X  -  6.3 

Y  ^  5.5 

CALL  P0RIG<X7Y) 
C 

C —    Scale  ay.es, 
C 

CALL  PSCftLE  (X_riATAr  lOJ-lr  1  fl .  0  r  FlRSTXr  nFI.  TAXy  rUVLNX) 

CALL  PSCALE  (Y^.DATA7  1021.  11.0.  FIRSTY7  DEI.TAY,  DiyLNY) 
C 

C —    Draw  X  axis. 
C 

CALL  PSLPEN(l) 

CALL  PAXIS(0.0.  O.Or  'R.-man  Sh.ift'r  '  rin  I  t- U.  I  '  r  220.  7?(n< 

-  ie.7  0.7  FIRSTX7  BELTAX.  nU'LNX.) 
C 

C —    Drsw  Y  axir.. 
C 

IF  (PLOT-NUMBER. NF.1)  THFN 

STR    -     '  Aitipl  i  t,Mdc  ' 
ELSE 

STR    ~     'Dotectar    Beei'^iori' 
END    IF 
CALL    PAXIS(0.0.     0.0.    RTR,     '     ',     120.     1201.     11..     V0.» 

-  FIRSTY.  DELTAYj  DiyLNY) 
C 

C —    Plot  data. 
C 

FIRDELd  )  -  FIRSTX 
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FIRDEl.  (3)  ^  DLL  TAX 
FIRriEL.(3)  -  FIRSTY 
FIRtiFL  <1)     ^    IiFLTAY 

CALL    F'SLPEN(2) 

CALL  PLINF  (X..tiATAF  Y.nATAr  10.?'lr  FIRPFL;  Or  '  '  ;  niVLNXr  tilVLMY) 

Plot  the  r.  i  a  n  2 1  - 1  o  -  n  o  i  G  c  r  i3 1  i  o  r.  n  t  t  h  p  i  r  r  p  f;  r  p  c  t  i  v  p  r-  u  1  r;  e  r  n  r-  i  t  i  o  n  . 

IF  (PLOT.NUMRFR.NF. 1 )  THFN 
DO  T  ^  1  )  15 

IF  (PNR(I) .NF.O. )  THFN 

Calculation  of  h o  r  i  r  o n t  a  .1  r  o n-  i  t  i  o 1 1 . 

X  ^  (MFANd  J/llJOO.  )«18. 

Select  vertical  porition  .:!rrorci'.i  n-i  to  uhpthrr  rulr,c  numhcr 
ir>  oven  or  odd. 

IF  (n«(I/?) .NF.I)  THEN 

Y  ^  11.7 
ELSE 

Y  ^  11.2 
ENDIF 


CALL  PPI  DTCXr  Yr  0) 

CALL  PFLTMM(SNR(I) >  S,  2) 

Show    ror.it ion    of    oridinal     r-ulnr. 

Y    ^    1  .0 

CALL  PPLOTCXt  Y.  0) 
CilLL  PFIXNM<l7  2) 
CALL  PPLOT(X.  0.0?  0) 
CALL  PPLOTR(0.0>  3.0,  1) 
END  IF 
END  DO 
ENDIF 

Wait  for  csrricae  rrturn  to  rlcDt-  the  nrrccn. 

CALL  BETUTXd,  '  '.  1.  INARY-  IBOT) 
CALL  FCLR5C 
CALL  FCLOSP 


RETURN 
END 
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c 

C        .SUBROUTINE  FFT 

C 

C        MAX-ll  FORTRAN  SOURCE  FILFNrtMF:  FFT.FOR 

c 

C  DEF'ARTMFNT    OF    FFCF  KYiNBrtS    STATE    UNIVI.-RSITY 

C 

C  REFERENCE!       Ni3<5ir    fthmcd    nnd    T.     N.Tt..-?r,-;,i.i!ri  t     Hi  Gcrpf.p-Timc    Sjlnclr.    rsrirj 

C  S<<r.t,o»r,j    Roston    PuhlishinS    Commny;     t9f)3.       pp.     tAO-161. 

C 

C***********«««*««««**««*«*«i(e*«*)!(**«)!!«*««««**«««S»««)!t«i![»)|(««««« 

C 

C        CALLING  SEQUENCE!   CALL  FFT(X»  Ni  INV) 

C 

C  PURPOSE!       This    routinp    performs    .t    (drw^rd    or    invcrr-r    ticci  m.^t  ion -in- 

C  f  roQuoncM    FFT    on    the    <3!?auencf'    X. 

C 

C        ROUTINE(S)  CALLED  OR  ACCFSSFH  PY  THIS  ROUTINE! 

C  None. 

C 

C        ARGUMENTS  REOUIRED! 

C 

C  X  r  o  m  r  .1  o  X 

C  The    pvrn<i    c o n t. b i  n  i  n •<    t h p    ?; g ct u p n r o    on    w h  i  r ^l    thr    FFT 

C  i  r>     to    hf?    riDrforiiiQd )     .TTr>o    tlir?    tirrny    in    uihjrh    thr'     result 

C  i  s    r  c  t  u  r  n  o  d  . 

C 

C  N  intcrtcr 

C  T  h  t?    1  o  n  a  t,  h    o  f    t  h  o    .3  r  r  3 '.)    X  . 

C 

C  INy  intcdcr 

-  Indicates    w h p t h p r    a    forward    (  0  )    or    i n v p r  r  p    (  .1  )     t  r n n r, - 

C  formir,    t oho    r'5?r formed. 

C 

SUBROUTINE  FFTtX,  H,     INW) 
COMPLEX  X(N) I  U;  T 
ITER  ^    0 
I  REM  -  N 
10     IREM  -  IREM/2 

IF  (IREH.EO.O)  GO  TO  20 
ITER  -:    ITER+1 
60  TO  10 
20     CONTINUE 
S  =  -1 

IF  (INV.EO.l)  S-=l 
NXPS-^N 

DO  50  IT-lr  ITER 
NXP-NXP2 
NXP2-NXP/2 

WPUR  -^  3.111592/FL0AT(NXP2) 
DO  10/1^1,  NXP2 

ARG  -  FLOAT <«-1 ) «UPUR 
U  -  CMPLX<C0S(^RH1  )  S*SIN(ARi5)) 
DO  10  MXP  -  NXP.  Nt  NXP 
Jl  -  MXP-NXP+M 
J2  --  J1  +  NXP2 
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X(  Jl  )-X(  11 


55 
60 


65 


70 
75 


X< Jl)+X( J2) 


GO  TO  55 


X(  Jl) 
40  X(J2) 

50     CONTINUE 
N2-N/2 
N1=N-1 
J  =  l 
no  65  I=-l.  Nl 

IF  (I. HE. J) 

T  -  X(J) 

X  <  J )  -  X  <  I ) 

X(I)  -  T 

K  ^-    N2 

IF  (K.GE.J)  GD  TO  65 

J  -  J-K 

K  =  K/2 

GO  TO  60 
J  ^  J  +  K 

IF  <INV.FQ.l)  GO  TO  75 
no  70  T  =  1  ,  N 
X(I)  =•  X(I)/FI.  OAT(N) 
CONTINUE 
RETURN 


END 
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The  objective  was  to  develop  a  method  for  the  automatic 
detection  of  signals  resembling  those  encountered  in  Raman 
spectroscopy.  First  a  matched  filter  is  used  to  enhance  the, 
spectral  peaks.  This  is  followed  by  a  detector,  which  was 
developed  through  the  application  of  statistical  decision 
theory. 

Decision  theory  involves  the  use  of  some  decision 
criterion  to  decide  between  two  or  more  alternative 
hypotheses.  The  Neyman-Pearson  criterion,  and  a  modifica- 
tion known  as  Wald  or  sequential  testing,  were  chosen  for 
study  because  they  require  relatively  little  knowledge  of 
the  signal;  in  particular  the  a  priori  probabilities  for  the 
signal  do  not  need  to  be  known.  The  Wald  test,  which  uses  a 
variable  number  of  samples  to  make  each  decision,  was 
finally  chosen  for  experimentation. 

For  the  matched  filter  the  impulse  respone  matches  that 
of  a  Lorentzian  peak  of  typical  width.  In  addition  to  the 
matched  filter  and  detector,  it  was  found  that  it  was 
necessary  to  fit  and  remove  a  trend  from  the  output  of  the 
filter  if  the  threshold  comparison  procedure  of  the  detector 
was  to  work  best. 

A  program  was  written  to  generate  signal  sequences  and 
perform  the  filtering  and  the  detection  procedure.  In  the 
detector  there  are  some  parameters  that  need  to  be  chosen  by 
the  observer.  These  are  the  theoretical  probabilities  of 
false  detection  and  false  dismissal,  and  also  what  is  termed 
the  assumed  signal  level.   Tests  were  made  on  several 
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sequences  of  signals  to  determine  how  well  the  system 
detects  signals,  and  also  to  determine  what  effects  the 
detector  parameters  have.  The  results  from  tests  on  four 
sequences  are  presented.  For  some  of  the  sequences  several 
tests  were  made  with  different  sets  of  parameter  values. 

The  results  show  that  it  is  possible  to  detect  peaks 
that  are  as  low  as  0.5  or  lower  in  signal-to-noise  ratio. 
However,  in  some  cases  peaks  with  signal-to-noise  ratios  as 
high  as  1.5  to  2.0  were  impossible  to  detect,  mainly  because 
of  the  effect  of  the  trend  removal.  Detection  at  low 
signal-to-noise  ratios  requires  accepting  an  increased 
likelihood  of  false  detections.  This  is  affected  primarily 
by  the  choice  of  the  theoretical  false-detection 
probability;  the  choice  of  the  false-dismissal  probability 
and  the  assumed  signal  level  did  not  appear  to  have  a 
consistent  effect  on  the  actual  probability  of  detection  or 
the  probability  of  false  detection. 

In  conclusion  it  is  suggested  that  the  system  may  be  a 
useful  addition  to  the  methods  of  peak  enhancement  that  are 
now  employed,  although  it  has  some  limitations  in  that  an 
estimate  of  the  noise  statistics  must  be  available,  and  also 
the  performance  may  depend  on  the  particular  form  of  the 
spectra. 


